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We propose that the solar neutrino deficit may be due to oscillations of mass-varying neutrinos 
(MaVaNs). This scenario elucidates solar neutrino data beautifully while remaining comfortably 
compatible with atmospheric neutrino and K2K data and with reactor antineutrino data at short 
and long baselines (from CHOOZ and KamLAND). We find that the survival probability of solar 
MaVaNs is independent of how the suppression of neutrino mass caused by the acceleron-matter 
couplings varies with density. Measurements of MeV and lower energy solar neutrinos will provide 
a rigorous test of the idea. 
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Introduction. Neutrino oscillation experiments have 
conclusively demonstrated that neutrinos have mass. 
Also, evidence has mounted that the expansion of our 
universe is in an accelerating phase caused by a negative 
pressure component called dark energy. While these two 
seemingly disparate advances are unequivocally among 
the most important of the last few years, our knowledge 
of both is woefully incomplete. 

Solar neutrino data provide the only evidence of matter 
effects on neutrino oscillations. The Large Mixing Angle 
(LMA) solution, while favored by reactor antineutrino 
data is somewhat discrepant with solar data, thus calling 
into question how well neutrino-matter interactions are 
understood. It is crucial that our understanding of how 
neutrinos interact with matter be confirmed or modified. 

Dark energy is troubling because the acceleration of 
the universe is a very recent phenomenon in its expan- 
sion history. This "cosmic coincidence" problem can be 
expressed as follows: Why are the dark matter and dark 
energy densities comparable today even though their ra- 
tio scales as ~ 1/a 3 (where a is the scale factor)? 

The coincidence that the scale of dark energy (2 x 10 ~ 3 
eV) 4 is similar to the scale of neutrino mass-squared dif- 
ferences (0.01 eV) 2 was exploited recently in Refs. 0, 3] 
to solve the coincidence problem. The authors of Ref. 2j 
considered the possibility of coupling neutrinos to dark 
energy by supposing that the dark energy density is a 
function of neutrino mass and imposing the condition 
that the total energy density of neutrinos and dark en- 
ergy remain stationary under variations in neutrino mass. 
Then neutrino masses vary in such a way that the neu- 
trino energy density and the dark energy density are re- 
lated over a wide range of a. 

A way to make the dark energy density neutrino-mass- 
dependent is to introduce a Yukawa coupling between 
a sterile neutrino and a light scalar field (similar to 
quintessence) called the acceleron. At energy scales be- 
low the sterile neutrino mass, the effective potential of 
the acceleron at late times receives a contribution equal 
to m v n v , where m v and n v are the active neutrino mass 
and number density, respectively. Supersymmetric mod- 
els of neutrino dark energy have been constructed Q. 



Model-independent tests of neutrino dark energy are 
cosmological 00. A strict relationship between the dark 
energy equation of state and neutrino mass is predicted. 
Further, since neutrino masses are predicted to scale with 
redshift approximately as (l + z)~ 3 in the nonrelativistic 
regime, cosmological and terrestrial probes of neutrino 
mass could give conflicting results. If the acceleron cou- 
ples both to neutrinos and matter, it may be possible to 
investigate this scenario through neutrino oscillations |]| . 
The coupling to matter is model-dependent. The effec- 
tive neutrino mass in matter is altered by the interactions 
via the scalar which in turn modifies neutrino oscillations. 

For environments of approximately constant matter 
density, a satisfactory approach is to parameterize the ef- 
fects of the nonstandard interactions by effective masses 
and mixings in the medium pj- However, for solar neu- 
trino oscillations it is not possible to account for the ex- 
otic matter effects by a constant shift in the oscillation 
parameters because the matter density in the sun varies 
by several orders of magnitude. 

In this letter we investigate solar MaVaN oscillations; 
these have not been studied previously. We will show that 
since the neutrinos propagate adiabatically, the specific 
dependence of the evolving masses on the acceleron po- 
tential is irrelevant, so the predicted survival probability 
depends only on the masses at their production sites. We 
then demonstrate how MaVaNs improve the agreement 
with solar neutrino data while being perfectly consistent 
with KamLAND data [jj. Finally, we illustrate via a 
calculation of the survival probabilities of atmospheric 
muon neutrinos crossing the earth's core that the scheme 
is consistent with atmospheric neutrino data. Since we 
focus on astrophysical and terrestrial neutrinos, the de- 
pendence of the neutrino mass on redshift is not pertinent 
to our considerations. 

Effect of acceleron interactions on neutrino 
masses. At low redshifts, the contribution to the neu- 
trino mass caused by the interactions of the acceleron 
with electrons and neutrinos can be written as 
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where \ Ui (A e ) is the Yukawa coupling of the acceleron to 



Vi (the electron). Throughout, when we quote values for 
A, we mean |A|. In principle, the scalar <j> has a mass, m^, 
that depends on n e and the n Ui . This dependence is weak 
since the underlying assumption in obtaining Eq. Q is 
that 4> does not fluctuate signficantly from its background 
value in the current epoch. The number density of the 
cosmic neutrino background in one generation of neutri- 
nos and antincutrinos is n„ uB ~ 112 cm~ 3 ~ 10 -12 eV 3 , 
the number density of relativistic neutrinos in the back- 
ground frame is n r v f , and the electron number density is 
n e . Here, m Vi are neutrino masses in a background dom- 
inated environment. We assume the heaviest i/j to be 
0(0.05) eV in the present epoch, and that as a result of 
their nonnegligible velocities, the neutrino overdensity in 
the Milky Way from gravitational clustering can be ne- 
glected [3. Then, m Vi essentially represent the masses of 
terrestrial neutrinos in laboratory experiments like those 
measuring tritium beta decay. Since the neutrinos under 
consideration are light, we do not expect the instabilities 
of highly nonrelativistic neutrino dark energy || . 

In principle, we should include a nucleon-acceleron 
Yukawa coupling. Since the electron-acceleron and 
nucleon-acceleron couplings are arbitrary (within bounds 
from gravitational tests), we can parameterize their com- 
bined effect on M, through A e , although this is not rig- 
orously true for two reasons: (1) Conventional matter 
effects for active neutrino oscillations do not depend 
on the nucleon number density njv- (2) The n e and n m 
distributions in the sun do not have the same shape |10l | . 
Nonetheless, this simplification suffices for our purposes. 

Tests of the gravitational inverse square law require 
the coupling of a scalar to the square of the gluon field 
strength to be smaller than O.Olmjv /Mpi ~ 10~ 21 [llj |. 
where rriN is the nucleon mass. Since we have chosen to 
embody the effects of the couplings of the acceleron to 
the nucleons and electrons in A e , the latter bound applies 
to A e . In the region of the solar core where yp neutrinos 
are produced, n° e ~ GOA^/cm 3 ~ 10 11 eV 3 [lfl. (Here 
and henceforth, we denote the electron number density 
at the point of neutrino production by rig). Thus, for A e 
close to its upper bound, A e n® ~ 10~ 10 eV 3 . 

The cosmic neutrino background contributes negligibly 
to the mass shift even for A„ i of 0(1). The pp reaction 
creates neutrinos with the highest number density in the 
production region (~ 7 x 10 -8 eV 3 ) and lowest energies 
[E v ~ 0.3 MeV) of all other processes in the pp chain and 
CNO cycle. Thus, pp neutrinos have the highest possible 
m v n r v el I 'E v , which for m„ of 0(1) eV is at most n^" B . 
In sum, the dominant contribution to the mass shift at 
the creation point arises from the X e n e term. 

We require that some Mj be 0(1O -3 — 10 -2 ) eV at neu- 
trino production in the sun. Then, for an assumed 
of 0(lO~ n ) eV 2 , we need X Vi ~ 10~ 4 - 10~ 3 . For this 
range of X Vi , the cosmic neutrino contribution in Eq. Q 
is five to six orders of magnitude smaller than the elec- 
tron contribution, and the pp neutrino contribution is 



eight to nine orders of magnitude smaller. The cosmic 
neutrino background density becomes dominant only af- 
ter n e drops by about six orders of magnitude. This does 
not happen until neutrinos reach the surface of the sun. 
As the neutrinos leave the sun, m„ approaches its back- 
ground value. The choice \ Vi ~ 10~ 3 serves more than 
one purpose. In addition to fixing the maximum values of 



Mi, it ensures that 



can be neglected for the entire 



path of the neutrinos through the sun. 

Solar MaVaN oscillations. In the framework of 
the Standard Model (SM) with massive neutrinos and 
conventional neutrino-matter interactions, solar (atmo- 
spheric) neutrinos oscillate with <5m 2 ~ 8 x 10~ 5 eV 2 
and 9 S ~ tt/6 (|<5m 2 | - 0.002 eV 2 and 6 a ~ tt/4 0). 
In our notation, <5m 2 (<5m 2 ) is the solar (atmospheric) 
mass-squared difference and 9 S , 8 a and 9 X are the mixing 
angles conventionally denoted by #12, O23 and 6*13, re- 
spectively @. We also know that solar and atmospheric 
neutrino oscillations largely occur independently of each 
other because 9 X must be small from the nonobservance 
of v e — > Vfj, oscillations at the atmospheric scale. In fact, 
data from the CHOOZ experiment demand sin 2 6 X < 0.05 
at the 2(j C. L. in the conventional picture. 

With the additional freedom that the M, provide, there 
is no reason to believe that the three neutrino oscilla- 
tion dynamics factorizes into the dynamics of two two- 
neutrino subsystems. Nevertheless, since our purpose 
here is to show that MaVaN oscillations are consistent 
with solar and atmospheric neutrino data while obeying 
the CHOOZ bound, we are entitled to accomplish our 
goal via construction. A simplifying assumption is that 
the decoupling of solar and atmospheric neutrino oscilla- 
tions continues to hold for MaVaNs. Then, the CHOOZ 
bound is automatically satisfied and we need to demon- 
strate that the two neutrino framework is adequate for 
both neutrino anomalies. 

The evolution equations for solar MaVaN oscillations 
in the two-neutrino framework are 

(mx-MxW) 2 M 3 (r) 2 
M 3 (r) 2 (m 2 -M 2 (r)) 2 

Here, Mj is a linear combination of those in Eq. QJ, 
U is the usual 2x2 mixing matrix, E v is the neu- 
trino energy and A{r) = 2\f2GF°n e {r)E v = 1.52 x 
10 _7 eV 2 n e (r) E v (MeV) is the amplitude for v e — e for- 
ward scattering in matter with n e in units of JVa/chi 3 . 
For typical 8 B neutrinos {E v ~ 7 MeV) n° e ~ 100, for 7 Be 
neutrinos (E v ~ 0.9 MeV) n° e ~ 90, and for pp neutrinos 
(E v - 0.3 MeV) n° e ~ 60. The matter term A at the 
points of origin is about 10 -4 eV 2 , 10 -5 eV 2 and 10~ 6 
eV 2 for 8 B, 7 Be and pp neutrinos, respectively. With our 
choice of |M 2 | of 0(1O~ 5 - 10~ 4 ) eV 2 at neutrino pro- 
duction, we expect nonstandard matter effects to be of 
the same order as standard matter effects. 
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We adopt a matter dependence of the form, 

M^r) = fa(n e (r)/n° e ) k , (3) 

where k parameterizes the dependence of the neutrino 
mass on n e , and fa is the neutrino mass shift at the point 
of neutrino production. As noted above we expect k to 
be close to unity, but we shall show that a wider range of 
k is allowed by oscillation data. We have implicitly made 
the approximation that all neutrinos are created with the 
same values of fa irrespective of where in the sun they are 
produced. Since almost all solar neutrinos are produced 
within r < 0.2r Q , for which n e falls by about a factor of 
3 from its value at the center of the sun, we consider the 
approximation to be reasonable. 

We make the parameter choices fa = mi = 0, fa> = 
0.0077 eV, fa = z0.0022 eV, m 2 = 0.0089 eV and 
9 = 0.62. The value of 5m 2 in a background dominated 
environment is m| = 7.9 x 10~ 5 eV 2 . We will sometimes 
refer to fa as MaVaN parameters and as background 
parameters. As we show, this set of parameters is consis- 
tent with KamLAND data and improves the agreement 
with solar data. 

The evolution of the mass eigenstates as they travel 
through the sun is governed by 
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where A(r) is the magnitude of the mass-squared differ- 
ence of the eigenvalues of the matrix in square brackets 
in Eq. @ and 9 m is the effective mixing angle in matter. 
The value of 9 m at the creation point of the neutrino is 



cos 26»° 
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2fa\ sin 29 



B 2 + (m 2 - fa) 4 sin 2 29 



4l4(A° sm29 + f4) 

where B = A Q — (777,2 — fa) 2 cos29, which yields the 
standard result in the limit that fa, fa — + 0. With 
Q(r) = 4E , the condition for adiabatic evolu- 

tion [13 is Q > 1. 

In Fig. ^ we show how 9 m and Q~ l depend on r/r© 
for E u = 0.1,0.74,5 MeV with k = 1. We do not show 
the evolution of A since it is smooth throughout. Notice 
the step in 9 m at r/r© - 0.0035 for E v = 0.74 MeV. (The 
energy at which this occurs depends on the background 
and MaVaN parameters chosen). The step manifests it- 
self as a large spike in Q" 1 ; Q < 10 only in a 1 keV spread 
around 0.74 MeV. While adiabaticity is violently violated 
in this narrow range of energy, it is undetectable because 
experimental resolutions are much larger than 1 keV. 

For all practical purposes, the evolution is adiabatic 
and the survival probability is given by the standard for- 
mula ,13| , P(v e -> v e ) = (1 + cos26>°„cos26>)/2, with 
008 26^ from Eq. Thus, we find that the survival 

probability of solar neutrinos is independent of k so long 
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FIG. 1: (a) 8m and (b) Q 1 as a function of t/tq for three 
representative energies. The adiabatic condition Q 2> 1 is 
violated for E v = 0.74 MeV at r/r Q ~ 0.0035 because d8 m /dr 
becomes very large. 



as the neutrinos propagate adiabatically. The depen- 
dence on the acceleron-matter couplings enters only at 
the production point of the neutrino via the fa. 

MaVaN oscillations vs. data. We now compare 
the predictions of this framework with solar data. To 
this end, we use the recently extracted average survival 
probabilities of the low energy (pp), intermediate energy 
( 7 Be, pep, 15 0, and 13 N) and hi gh energy ( 8 B and hep) 
neutrinos; for details see Ref. PJ- From Fig. we 
see that the MaVaN survival probability almost passes 
through the central values of the three data points. 
The agreement with intermediate energy data is remark- 
ably improved compared to the LMA solution because 
PMaVaN^e v e ) approaches sin 2 9 for lower E v than 
for PsM^ye — * v e ). For the same solution, it is possible 
for pp neutrinos to have a higher survival probability than 
the vacuum value, 1 — 0.5 sin 2 29. This control over the 
width of the transition region and the larger difference 
between the survival probabilities of the pp and 8 B neu- 
trinos (than cos 2 9 cos 29 for the LMA solution) is a result 
of the freedom provided by the additional free parameter 
''fa. Keeping in mind that the survival probability of the 
neutrinos incident on earth is independent of A;, if we set 
fa — k — 0, we recover the standard MSW case with 777,2 
replaced by 7772 — fa. The fc-dependence reappears for 
neutrinos passing through the earth. 





0.8 




0.6 


T 




f? 


0.4 


a. 






0.2 









10 





"---JJvIA 









10 u 

Ey [MeV] 



10' 



FIG. 



2: P(v e —* u e ) vs. E v for MaVaN oscillations (solid 
curve). The dashed curve corresponds to conventional oscil- 
lations with the best-fit solution to KamLAND data. 

An important question is whether MaVaN oscillations 
are consistent with KamLAND data. The solid curve in 
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FIG. 3: (a) dP(v e — > i/ e ) for reactor antineutrinos incident at 
KamLAND (solid), and dP{v2 — > z^ e ) for solar neutrinos pass- 
ing through the center of the earth (dashed), (b) dPiv^ — » f M ) 
for atmospheric neutrinos passing through the earth's core. 



Fig.Et is dP{v e — > V e ) = PsM{Ve — * V e )- P M aVaN{Ve ~ > 

P e ), for a mean KamLAND baseline of 180 km and energy 
resolution 7.3%/ \J E(AIeV). Here, Psm is calculated for 
k = 1/2 with 6m 2 = 8 x 10~ 5 eV 2 and = 0.55; the 
latter are the vacuum parameters favored by KamLAND 
data. Since mass- varying effects in the earth scale like the 
ratio of electron number density in the earth to that in 
the sun, the effects are larger for smaller k. Any k > 1/2 
produces a dP that is acceptable. 

Another relevant question is if earth-matter effects are 
substantial for solar MaVaNs. Since v e with energy above 
a few MeV exit the sun and arrive at the earth in the 
intermediate neutrino mass eigenstate (denoted by ^2), 
it is appropriate to study V2 —> v e transitions to assess the 
size of these effects. The dashed curve in Fig. OK is the 
energy-averaged dP{i>i — > v e ) (assuming a 10% energy 
resolution) for neutrinos passing through the center of 
the earth in which case matter-effects are expected to be 
enhanced. PuaVaN(^2 ~ * v e) deviates only slightly from 
the usual matter oscillations and is in accord with a tiny 
day-night effect as required by Super-Kamiokande and 
SNO data. 

We next show in Fig. that atmospheric neutrino 
data are also consistent with 2-neutrino — ► v T oscilla- 
tions for MaVaN parameters of comparable size to those 
in the solar sector (k = 1/2, /if = 0.01 cV, = 0.003 
eV and 7713 = 0.047 eV). Here, Psm is calculated for 
8m 2 = 0.0021 eV 2 and 9 = tt/4. dP{v^ -> VfJ ) is aver- 
aged over the earth's core (cosOz = 0.8 — 1, where 0z is 
the nadir angle), for a 10% energy resolution. We also 
confirm that dP(v^ —> i/„) at the K2K baseline is well 
below experimental sensitivity. 

Conclusions. We have shown that oscillations of vari- 
able mass neutrinos (that result in exotic matter effects 
of the same size as standard matter effects) lead to an im- 
proved agreement (relative to conventional oscillations) 
with solar neutrino data while remaining compatible with 
KamLAND, CHOOZ, K2K and atmospheric data. 

MaVaN oscillations are perfectly compatible with solar 
data because the survival probability can change from 
a higher-than- vacuum value (at low energies) to sin 2 9 
(at high energies) over a very narrow range of energies. 
Since the neutrino propagation is highly adiabatic, the 



survival probability of solar neutrinos is independent of 
k as defined in Eq. 

Whether or not an explanation of solar neutrino data 
requires MaVaN oscillations will be answered by experi- 
ments that will measure the survival probability of MeV 
and lower energy neutrinos. As shown in Ref. [15j, other 
tests in reactor and long-baseline experiments emerge 
when the scheme presented here is embedded in a com- 
prehensive model that can explain all extant neutrino os- 
cillation data including the LSND anomaly and a future 
MiniBooNE result. 

We have considered neutrinos with background mass 
of 0(0.01) eV. For such light neutrinos, only model- 
dependent (neutrino oscillation) tests of the MaVaN sce- 
nario are viable because the model-independent (cosmo- 
logical) tests become inoperable. There are two reasons 
for this: (1) The dark energy behaves almost exactly as 
a cosmological constant today. (2) If these light neutri- 
nos do not cluster sufficiently, the local neutrino mass 
is the same as the background value, which is below 
the sensitivity of tritium beta-decay experiments. Then, 
high-redshift cosmological data (which should show no 
evidence for neutrino mass) and data from tritium beta- 
decay experiments will be consistent. 
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